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ABSTRACT 

The dark cloud Lynds 1622 is one of a few specific sites in the Galaxy where, 
relative to observed free-free and vibrational dust emission, there is a clear excess 
of microwave emission. In order to constrain models for this microwave emission, 
and to better establish the contribution which it might make to ongoing and near- 
future microwave background polarization experiments, we have used the Green 
Bank Telescope to search for linear polarization at 9.65 Ghz towards Lynds 1622. 
We place a 95.4% upper limit of 88 fiK (123 fiK at 99.7% confidence) on the total 
linear polarization of this source averaged over a 1'.3 FWHM beam. Relative to 
the observed level of anomalous emission in Stokes I these limits correspond to 
fractional linear polarizations of 2.7% and 3.5%. 

Subject headings: ISM: dust, ISM: individual (Lynds 1622), polarization, radio 
continuum: ISM 
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Introduction 



The discovery of anomalous dust-c orrelated emission 
Background exper iments in the 1990s (ILeitch et al 



Kogut et al. 



19961 ) led to a recalculation by 



1997 



jy multi-frequency Microway e 



de O 



Draine & Lazarian 



iveira-Costa et al. 



1998a 



dipole emission from small rotating dust grains, first considered by 



19971 : 



) of the electric 



Ericksonl (119571 ). The 



emission is spatially corr elated with dust, but 



of the thermal emission (IFinkbeiner et al. 



emission, or "Foreground X" (jde Oliveira-Costa et al. 



ar in excess of a reasonable extrapolation 



1999) and is therefore called "anomalous" dust 



2002) becaus e of its interference 



Draine fe Lazarian 



with co smic microwave background (CMB) experiments. Since the 
( 11998a! ) paper on spinning dust, and a later paper proposing magnetic dipole emission 



from Fe containin g grains ( 



been refined (e.g. 



Draine fe Lazarian 



de Oliveira-Costa et al. 



1998 



1999) measurements of t his component have 



Finkbeiner et al. 



2002 



Bandav et al. 



2003 



de Oliveira-Costa et al. 



2004 



Finkbeiner et al 



2004 



Davies et al. 



20061 ) using various data 



sets. 

One yet untested prediction of the spinning dust theory is that the emission should 
be moderately polarized (3 — 7%) at low frequencies, falling to under 2% above 20 GHz 



(ILazarian fe Draine 



2000j). At 10 GHz 3-5% linear polarization is expected. In contrast 



magnetic dipole fluctuations in larger single-domain grains would sh ow a strong linear 



polar ization (~ 10%) below 10 GHz, rising to over 30% at 100 GHz ( jDraine fe Lazarian 
19991 ). although the details depend on the shape, iron fraction, and magnetic domain 
configuration of the grains. Because the polarization properties of the anomalous dust 
emission are unknown and are of immense interest to the CMB community, a deep exposure 
on a cloud known to contain the anomalous emission would make a substantial contribution 
to our understanding of this component. 



Although anomalous, dust-correlated microwave emission has been detected 
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statistically by numerous experiments, at the time of writing the number of indivi dual 



2005|). 



Finkbeiner et al 



Battistelli et al 



2002; 



Casassus et al 



lines of sight on which it is seen a r e but three: the NC P loop (ILeitch et al 
Lynd s 1622 ( 



20041 ); and Perseus ( IWatson et al 



1997) 



(120061 ) report the only measurement of the polarization of the 
anomalous emission, a detection of 3% fractional polarization in the Perseus cloud with the 
COSMOSOMAS experiment. We have used the GBT at 9 GHz to search for polarization 



towards Lynds 1622 (L1622). 



2. Observations & Calibration 

We chose to search for polarized emission from L1622 along the line of sight 05:54:23, 
+01:46:54 (J2000); this is the same pointing position used by Finkbeiner et al. in previous 



9 GHz measurements of LI 622. 



observed by CBI (jCasassus et al. 



t is also near the peak of 30 GHz excess emission 



20061 ). We observe at 9 GHz, a frequency where the 
product of fractional linear polarization and total intensity for spinning dust is expected 
to be detectable (a few mK). This is also the lowest frequency GBT receiver equipped 
with circularly polarized feeds, which are essential to obtaining the broadband continuum 
stability needed for the Stokes Q and U measurements. Two sub-bands were chosen by 
examining site RFI monitor data, each of 200 MHz total bandwidth, centered on 8.65 and 
9.65 GHz, respectively. Much wider bandwidth measurements are in principle possible, but 
RFI was a key consideration. The receiver temperature is smooth and low in these regions 
and avoid known spectral resonances. 

The GBT X-band receiver has a single, dual-circularly polarized feed horn. The GBT 
Spectrometer was used to form all four auto- and cross-correlations between the IF signals 
from each of LCP and RCP, providing in principle instantaneous measurements of all four 
Stokes parameters (/, Q, U, V). The Stokes parameters describing linear polarization (Q, U) 
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principally comprise combinations of the cross- correlations, which has the advantage that 
receiver gain fluctuations and atmospheric emission variations are suppressed. 

Successive pairs of On/Off measurements on L1622 were performed, each lasting 48 
seconds and consisting of 1-second integrations. The On/Off pairs were matched in hour 
angle to provide approximately the same track in azimuth and elevation for each pair; 
allowing for scan-related observing overheads we found we required a separation of ~ 66 
seconds in Right Ascension. This strategy minimizes the potential effect of polarized ground 
spillover, which would be expected to have a signature that changes with azimuth and 
elevation. Note however that the expected level of spillover is low due to the clear aperture, 
off-axis design of the GBT; deep integrations with the GBT 26 — 40 GHz receiver show that 
the level of any systematic signal in total intensity is < 100 /xJy. The choice of trailing field 
was constrained by a compromise between observing efficiency and the benefits of short 
cycle times, particularly for the total intensity measurement. The timing accuracy of the 
GBT control system was not sufficient to allow an analogous LEAD region. In all 700 
On/Off pairs were collected in 30 hours of observing. Every hour the telescope pointing and 
focus is checked on the nearby calibrator 3C138. 3C138 is also a well-measured polarization 
calibrator and we perform an On/Off measurement of 3C138 with the system configured 
identically to our L1622 polarization observations. 

Prior to any calibration the data are flagged for RFI. Several fixed frequency ranges 
show common interfence events and are flagged in all scans (9.57579 — 9.5777 GHz; 
9.6673 - 9.6688 GHz; 9.7005 - 9.7018 GHz; and any data over 9.730 GHz). Beyond this the 
RL and LR data for each IF are searched and integration with a > 5cr spectral feature is 
flagged for rejection. 



The ful 



Heilesl fcOQlh 



-Stokes data were ca 



Heiles et al. 



fl2001b 



ibrated using a variation on the procedures described in 



aj), which we briefly review. A single noise diode is used 
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to inject noise into each of LCP and RCP, and this (coherent) signal is used to determine 
the phase difference between the two IF paths. The dominant contributor to the R/L 
phase is the path length difference between the signal chains (~ 70 cm). The response 
of the instrument to celestial polarization is described by the Mueller matrix (MM). We 
determine the MM from the 3C138 measurements assuming a fractional linear polarization 
of 11% and a parallactic angle of 170° of the linear polarization pseudovector. The data are 
also amplitude-calibrated using 3C138, for which we assume flux densities 5*8.65 = 2.43 Jy 
and 5*9.65 = 2.20 Jy. This allows the determination of an effective Stokes I flux density of 
the calibration diode, S ca i, as a function of frequency across our observing band. After 
correcting for the measured LR phase difference and calibrating data for the bandpass as a 
function of frequency from the S ca i data, means across the band are taken to form nominal 
Stokes (J, Q, U, V) measurements. Mueller Matrices are computed from these data on the 
calibrator, and applied to the equivalent data on L1622. The result of these procedures 
is a position-switched Stokes (/, Q, U, V) in Janskys per beam for each of 720 ON-OFF 
measurements of L1622 in our dataset. This is converted into a main-beam filling equivalent 
surface brightness by multiplying by the GBT gain (taken to be 1.95 K/Jy at 9.65 GHz) 
and dividing by the main beam efficiency (r] B = Q mb /Q ant = 93%). 

The beam properties were determined by scans along six evenly spaced directions 
through 3C286, and orthogonal (az/el) scans through 3C138. We adopt a beam width of 
1'.30 (FWHM). 

The noise level was estimated from the data for each measurement by a robust, iterative 
procedure. First the the median absolute deviation^ of the measurements is computed for a 

x The median absolute deviation (MAD) is defined as Median(\xi — Median(xi)\) and is 
much less sensitive to outliers than the variance, although its sampling variance is greater. 
For a Gaussian distribution MAD = 0.674a; in data analysis, MAD values are normalized 
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Fig. 1. — Measurements of 3C138 (black pluses) and the fitted Mueller Matrix describing the 
GBT polarization response (green triangles). The quantities shown are, from top left, the 
average autocorrelation (LL+RR); the real part of the LR cross correlation; the imaginary 
part of the LR cross correlation; the the difference between the autocorrelations. 

given IF and Stokes parameter in a 3- hour buffer centered on the time of the measurement. 
We refer to the scatter so computed in this first pass on the data as <j\. Data more than 
6o"i from the mean are rejected iteratively, i.e., if a data point (s) over the threshold exists, 
the worst outlier is rejected, the a recomputed, and the process repeated until there are 
no data over the threshold. This rejects ~ 2% of the data. After this process is completed 



to a Gaussian-equivalent a. 
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the noises are recomputed using the variance in a sliding 3-hour buffer, resulting in a 
noise estimate Oi on the outlier-rejected data. Data with cr 2 values more than 5 times the 
expected radiometer noise (S maX: Qu = 5 x 0.2 mJy) in Q and U are rejected. Stokes / 
data with noises more than 5 times the best RMS (S max j = 5x10 mJy) are rejected. This 
removes data collected in periods of less than optimal weather. Final results are computed 
as the weighted mean of individual measurements, with weights of o%T 2 . The final results for 
each polarization are shown in Table [2] along with the average noise levels. Results for the 
8.65 GHz channel are consistent with those from the 9.65 GHz channel but the Stokes Q 
and U parameters show noise levels a factor of ~ 2 higher, perhaps reflecting a hardware 
instability in these channels; we exclude them from our final results. The large variation in 
noise levels between the Stokes J, Q, U, and V results is explained by the fact that with 
our observing technique both Stokes I and V are affected by receiver gain fluctuations, and 
Stokes I is additionally affected by fluctuations in atmospheric emission, whereas Stokes Q 
and U largely are not. 

A histogram of the individual measurements, normalized by the oi value for each 
measurement is shown in Figure [2] along with the best-fit Gaussian to the distribution which 
should have a Gaussian a close to unity if our noise estimate is accurate. For Stokes Q and 
U the fitted Gaussians are within 5% of unity; for / and V, within 12% of unity. While the 
data on the whole are well-described by a Gaussian distribution there are a small number 
of outliers. To determine the sensitivity of our results to these outliers, and to the noise 
estimate that results from our pipeline, we have varied our data filtering parameters in a 
suite of tests summarized in Table EJ The controlling parameters of our filtering procedure 
were systematically varied; the fully-averaged 9.65 GHz Stokes parameters that result for 
each variation are shown in Figure [3j Note that our adopted parameters are denoted as 
"Test 0". 
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The data binned by time and Parallactic Angle (PA) are shown in Figures H] and 
respectively. The overall consistency of the data are good, with the most significant 
deviation in Stokes V. The absence of significant variations in Stokes I, Q, and U with PA 
indicates that any residual ground signal is lower than the sensitivity achieved. 





Fig. 2. — The distribution of individual measurements divided by their noise estimates (8.65 
GHz is black and 9.65 GHz is magenta). The dashed lines show the best-fitting Gaussian 
to each distribution. In all cases the Gaussian has a a within 4% of unity indicating the 
accuracy of our noise estimate. 

Joint la, 2a, and 3a confidence intervals on the 9.65 GHz linear polarization signal are 
shown in Figure EJ These represent the difference in Stokes Q and U between the ON and 



-10- 



10 
5 


-5 
-10 



10 



1 5 



20 



25 



Test # 



o -0 



05 
00 
05 



4 + ) I f ^ ~ - ^ + + + ' + + + f + I 4 4- 4- - • • r i + 



15 



20 



25 



-0 



+ + + 

10 



Test # 



+ + + 



Test § 



I + ' 

15 



-±-i-i-i-i"i-*-i"f- - - 



20 



25 



2 












* 1 














+ + + + , _ - - - 


" + 1 + + + 


+ + | + + + + + 




> -1 










-2 












5 


10 


15 20 


25 





Test § 



Fig. 3. — Fully averaged Stokes I, Q, U, and V results for a range of filtering parameters. 
The dashed line shows our adopted value (test 0). 

OFF source locations, and on the assumption that the background signal is unpolarized at 
our observing frequency are limits on the polarization on the line of sight to L1622. 

We note that there is a negative signal in Stokes I on the line of sight to L1622. This 
is a robust result, likely due to small-scale variations in the free- free signal in the region; we 
discuss this possibility further in § I3.ll 
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Fig. 4. — Mean values for each of Stokes I,Q,U, and V with data binned in time over the 
course of the observing run. 

3. Results & Conclusions 

3.1. Total Intensity Signal 

We have measured a total intensity signal (On-Off) of — 5.1mK (— 2.6mJy/Bm) 
towards L1622, indicating that at 9.65 GHz the sky is brighter off-source than on-source. 
To understand the nature of this signal we must consider both spinning dust and free-free, 
which will each be significant at 9 GHz. Fortunately high-resolution Ha and 100 jum data 
that cover the region of interest are available. 
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Fig. 5. — Mean values for each of Stokes I,Q,U, and V with data binned in parallactic angle. 



L1622 is embedded within the Orion star forming comp 



ex, in a vicinity of significant 



20011 ) we determine that 



free-free emission. Using the SHASSA map ( iGaustad et al. 
there is a 20 Rayleigh difference between our on-source and off-source locations. The 
Ha measurements suffer from significant extinction by dust, which must be accounted 
for in or der to accurate l y pred ict the radio free-free signal. We estimate the extinction 



from the 



Schlegel et al. 



(119981 ) maps of temperature-corrected 100/im dust intensity. 



The SHASSA and IRAS 100 jttm maps are shown in Figure [7J along with our On and 



Off source positions. For dis p 



ay p urposes we use the more recent and higher- quality 



Miville-Deschenes fc Lagachd (120051 ) maps, however, the dust extinction relations are 
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Fig. 6. — Joint 68.3%, 95.4%, and 99.7% confidence regions for Stokes Q and U. 



based on the 



Schlegel et al. 



(jl998l ) D T maps. In the SHASSA map L1622 appears as 
a dark cloud, ostensibly between us and the majority of the Ha emission. The peak 
dust-corrected infra-red brightness on-source is ~ 160MJy/Sr; allowing for a characteristic 
background level of 60 MJy/Sr in the vicinity, the peak emission due to th e dark cloud itse lf 



is 10 MJy/Sr. With an Ha extinction relation of 0.0462 mag/(MJy/Sr) (IDickinson et al 



20031 ) we find an overall Ha extinction of 4.6 mag. At this high level of absorption the 
free-free predictions are not reliable. We can nevertheless estimate the characteristic level 
of free-free emission variations nearby, which we find to be ±5 mK on scales of our beam 
throw (66 seconds or 16'. 5) using the T e = 70007^ scaling of Ha intensity to radio surface 
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brightness (60.9 /iK/Rayleigh at 10 GHz — Dickinson et all (120031 )) 



With the limited frequency coverage of our two-point measurement we are unable to 
distinguish directly between free-free and " anomalous" exces s emiss ion. To do this we can 



we use the 8 and 10 GHz measurements of 



Finkbeiner et al. 



(120021 ). who, in linear scans 



across L1622 d etermine a peak dust- correlated excess signal of 4 mK. We can compare this 



to the 31 GHz 



Casassus et al 



(120061 ) detection of excess anomalous emission in L1622, 



which yield a dust emissivity of 24.1 ± 0.7 pK/ (MJy/Sr) 



GHz with a 37% CNM + 63 % WNM (IDraine fc Lazarian 



at 31 G Hz). Extr apolating to 9.5 



2006 



Finkbeiner et al 



1998bl ) dust SED (ICasassus et al. 



2002J) gives a dust emissivity of 54.6 yuK/(MJy/Sr), consistent with 



the value de Oliveira-Costa et al. (1999) give at 10 GHz, 50 /xK/(MJy/Sr). These scalings, 
together with the extinction-corrected levels of dust emission in the field, predict a 9.65 
GHz spinning dust signal of 3 mK. We adopt an estimated 9.65 GHz,excess dust-correlated 
signal level in Stokes I of 3.5 ± 0.5 mK. 



Finkbeiner et al. 



(120021 ) use the same 



We note that this measurement and that of 
central position on L1622 but different reference (off-source) positions; the chop throws are 
comparable (16' vs 12'). 



3.2. Polarization Signal 



With the GBT at 9 GHz we have determined that Q = 35±30/zK and U = 26±34/xK. 
Figure [6] shows the joint two-parameter 68.3%, 95.4%, and 99.7% confidence regions for 
fractional Q and U. We employ a simple Maximum Likelihood approach to set a limit on 
the total linear polarization that avoids the Rice bias. The likelihood of a linear polarization 
p = \JQ 2 + U 2 and polarization angle 6 given statistically independent data Q b s , U obs is 



L(p, 8\Q obs , U obs ) oc Exp 



(Q-Qobs) 2 , (U-U obs ) 



2a? 



+ 



2a? 7 

Uobs 



(1) 
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Fig. 7. — Ha map (left) and 100 /zm map (right). Our on and off positions are shown as 
green circles. 

where Q = pcos 28 and U = psin28. Marginalizing over the position angle of polarization 
we have 



L(p\Q obs , U obs ) oc / ddExp 



(Q — Qobs) 2 



(U ~ Uobsf 

2o* 



(2) 



The PDF is normalized to unity, and confidence intervals are determined by integrating the 
PDF. We find 95.4% and 99. 7& upper limits on the total linear polarization p of 88 fiK 
and 123 fiK, respectively. For a nominal Stokes I dust signal of 3.5 ± 0.5 mK (§ 13. II) . and 
assuming that the free-free signal is unpolarized, these limits correspond to fractional linear 
polarizations of 2.7% and 3.5%. 



These limits are consistent with the results of Battistelli et al. on Perseus, who find a 
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fractional linear polarization of 3.4^}g%. If electric dipole emission is primarily responsible 
for the observed microwave excess, a population of small dust grains is required, with 
tyipcal radii of a few nm. These grains can be efficiently aligned by paramagnetic resonance 
relaxation (Lazarian & Draine 2000), which would give rise to up to ~ 5% fractional linear 
polarization at 9 GHz. Both our result and that of Battistelli et al. are lower than this, 
but probably within theoretical uncertainties. In contrast, magnetic dipole fluctuations in 
larger, single-domain ferrous grains will show a fractional linear polarization ~ 10% below 
10 GHz, flipping orientation and rising to as much as 33% at 100 GHz (Draine Sz Lazarian 
1999). Many of the magnetic dipole models are excluded by our measurements. 



Conclusions 



We have presented a 9.65 GHz limit on linear polarization towards the dark cloud 
Lynds 1622, a well- measured locus of anomalous Galactic microwave emission; the total 
degree of linear polarization is < 88 fiK at 2a (< 123 /iK at 3a). Assuming a 3.5 mK 
stokes I spinning dust signal, consistent with independent m easurements of L1622 at 8 



through 32 GHz (jFinkbeiner et al. 



2002: 



Casassus et al. 



20061 ). and that the free- free signal 



is unpolarized, we have limits on the degree of linear polarization of 2.7% and 3.5%. These 
limits are consistent with the expected linear polarization of small rotating grains, and 
with the 3% fractional linear polarization measured by Battistelli et al. towards Perseus, 
but inconsistent with many models for the anomalous emission based on magnetic dipole 
fluctuations in ferrous grains. Such low levels of linear polarization would also be unusual 
for soft, dust-correlated synchrotron, which has been s uggested as a possib le origin for 



2003 . For L1622, 



the generally observed anomalous microwave emission (IBennett et al. 
however, soft synchrotron was already strongly ruled out by existing total intensity data. 
In addition we see that the total intensity signal on the line of sight to L1622 is fainter 
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than on the nearby comparison (Off-source) region. The negative signal we observe is 
consistent with the level of free-free fluctuations in the region: considering the expected 
+3.5 mK signal, the observed -5 mK signal indicates an 8.5 mK gradient between our 
ON and OFF positions, consistent with our estimate of ±5 mK variations over the region 
from Ha maps. This underscores the need for good frequency coverage in attempts to 
charactize the spinning dust foreground at low microwave frequencies. The relatively weak 
polarization signal observed in L1622 and Perseus bodes well for future CMB polarization 
experiments. It should be appreciated that while these concentrated loci of emission are 
important first steps on the road to understanding the anomalous microwave emission, 
the physical conditions are quite different from those in the diffuse ISM, where foreground 
properties are most important for future CMB expreiments. To study these regions more 
sensitive observations are necessary. 

The National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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Average ON-OFF Typical Per-Measurement Noise 

I -5.1±1.4mK 42 mK 

Q 35±30./iK 0.81 mK 

U 26±34/iK 0.93 mK 

V -0.1±0.3mK 10. mK 

Table 2: Final results of 9.65 GHz full-Stokes polarization measurements of L1622 and typical 
per-measurement noise level. 





P.A. 


MJD 






P.T.E. 


Xl P-T.E. 




0.44 


78% 


0.74 56% 


Q 


0.98 


42% 


1.97 10% 


u 


0.52 


72% 


0.54 70% 


V 


2.44 


4% 


1.92 10% 



Table 3: xt values for data binned in Parallactic Angle and time, with respect to the final 
fully averaged values, all with v = 4. 
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